Modelling the Berezinskii-Kosterlitz-Thouless Transition in the NiGa 2 S4 
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In the two-dimensional superfluidity, the proliferation of the vortices and the anti- vortices results 
in a new class of phase transition, Berezinskii-Kosterlitz-Thouless (BKT) transition. This class 
of the phase transitions is also anticipated in the two-dimensional magnetic systems. However, its 
existence in the real magnetic systems still remains mysterious. Here we propose a phenomenological 
model to illustrate that the novel spin-freezing transition recently uncovered in the NMR experiment 
on the NiGa2S4 compound is the BKT- type. The novel spin- freezing state observed in the NiGa2S4 
possesses the power-law decayed spin correlation. 

PACS numbers: 75.30.Gw, 75.40.Cx, 75.40.Mg 
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As the thermodynamic conditions in the environment 
change, for example the pressure or the temperature, 
matter transforms from one state to another. Steam lev- 
itates from the top of the hot coffee; ice melts in the soft 
drink. The phase transition is ubiquitous in our daily 
life. Investigating the phase transition is always the cen- 
tral subject in physics. Most of the phase transitions 
can be understood by the distinct physical properties be- 
tween the phases. For example, in the vapor-liquid tran- 
sition, the unit volumes per mole of the molecules are 
different from the vapor to the liquid. In the ferromag- 
netic transition, the spins orientate randomly at the high 
temperature side but start to align along the same direc- 
tion resulting in the net magnetization at the low tem- 
perature side. However, in the two-dimensional super- 
fluid, the Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion that happens when vortices and the anti-vortices 
proliferate does not separate the phases with distinct 
thermodynamic quantities [lj. Instead, the superfluid- 
ity correlation changes from the power-law behavior to 
the exponentially decayed one as we across the tran- 
sition from the lower temperature side. In the two- 
dimensional magnetic systems, the BKT transition is also 
anticipated in all easy-plane Heisenberg models. On the 
other hand, the power-law decayed spin correlation may 
play important role in the high transition temperature 
superconductor Q. Therefore, understanding the spin 
dynamics in the critical phase has become tremendously 
important. 

NiGa2S4 is originally synthesized intending to realize 
the spin liquid proposed by Anderson few decades ago@. 
It is the layered material that the spin-1 Ni 2+ ions form 
the ideal two-dimensional triangular network. As an anti- 
ferromagnetic insulator, NiGa2S4 exhibits no long-ranged 
magnetic ordering down to 350mK shown in the specific 
heat capacity, the magnetic susceptibility, and the neu- 
tron scattering experiments At 1.5K, the Edwards- 
Anderson order parameter, Q = l/iV^V < Si > 2 , that 
tells the spin moment, is measured at 0.61, where N is 
the total number of spins. This vastly reduced spin mo- 
ment indicates the presence of the strong quantum fluc- 



tuation that is highly favorable by the scenario of the 
spin liquid. However, in the recent experiments of the 
Ga nuclear magnetic resonance (Ga-NMR) [Bj], a tem- 
perature Tf is found around 10K below which the spin 
dynamics is slow down and the freezing behavior is ob- 
served. As approaching the Tf from above, both the 
nuclear spin-lattice relaxation rate 1/Ti and the nuclear 
spin-spin relaxation rate I/T2 diverge. Moreover, spins 
do not freeze immediately at Tf but persist fluctuating 
down to 2K found in the nuclear quadruple resonance 
measurement (NQR). Below 2K, the Ga-NQR spectrum 
becomes very broad and featureless, which implies the 
formation of the static inhomogeneous internal magnetic 
field. First, intuitively, the static internal magnetic field 
occurs when the spins freeze up completely. In this case, 
the Edwards- Anderson order parameter should be close 
to the quantum number of the spin angular momentum. 
Second, all thermodynamic quantities change smoothly 
at the transition temperature Tf, but the 1/Ti and the 
I/T2 diverge in the NMR signals. Therefore, the later 
NMR experiment apparently looks inconsistent with the 
previous measurements. In this Letter, we shall provide 
a consistent picture to compromise all the experimental 
results. Most importantly, we illustrate that the novel 
spin-freezing transition at T=Tf is the long-sought BKT 
transition in the two-dimensional magnetic systems. 

The way to compromise with all the experiments is to 
consider the state that contains both the freezing spins 
and the fluctuating ones. Then, both of them can be 
observed simultaneously in the experiments. Before ex- 
plaining further, let us begin by reviewing the spin con- 
figuration, depicted in Fig. (la), observed in the neutron 
scattering experiment. The correlation has the wave 
vector (1/6, 1/6,0) with the wavelength 2n/3a, where a 
names the lattice constant between two Ni 2+ ions. This 
wave vector simply means that along the ai direction 
the periodicity is of 6 sites and so is it along the a 2 , and 
along the b = (1, 1,0), the periodicity is of 3 sites. For 
convenience, we highlight the triangles in the Fig. (la) 
as the reference. With respect to the reference trian- 
gles, spins have the " l-in-2-out" (one spin points in and 
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FIG. 1: (Color online) (a)The spin structure observed by the 
neutron scattering in the Ref. |4[. There are two arrows on 
every site. The colored one is the spin orientation, and the 
thin black one is the local easy axis, (b) The ground state of 
the anti-ferromagnetic quantum Ising model on the triangular 
lattice. "-" indicates the local state and the "+" is the 
local " + state. "0" is the linear superposition of the 
and the states. 

two spins point out with respect to the triangles) or the 
"2-in-l-out" configurations. Inspired by the observation 
in the NMR experiment that the spin correlation starts 
to develop below the Curie- Weiss temperature 80K, we 
assume the existence of an easy axis on every spin site 
and its orientation is either parallel or anti-parallel to the 
current spin configuration. In order to manifest the anti- 
ferromagnetic nature, we assign the local +z axis to be 
the " all-in" or the " all-out" with respect to the reference 
triangles alternating over the whole lattice. An example 
of the orientation of the local axes is depicted as the black 
thin arrows in the Fig. (la). In this transformed coordi- 
nate, spins are either +1 (parallel to the +z direction) or 
-1 (anti-parallel to the +z direction). 

Considering the quantum fluctuation explicitly, we 
propose the following phenomenological Hamiltonian in 
the transformed coordinate 

H = j"£at**-KY,<rf, U) 

<ij> i 

where a z take the eigenvalues +1 and -1, the J is the pos- 
itive phenomenological coupling constant, and the K is 
the positive measure. Since the a x have the off-diagonal 
matrix elements that connect the +1 and the -1 states, 
the K measures the spin-flipping process between the 
+1 and the -1 states. For small K, the ground state 
of the Eq.(l) on the triangular lattice is given in the 
Fig. (lb). Using our coordinate transformation, we can 
map the spin configuration in Fig. (la) to the ground 
state of Eq.([T]). If the spin is parallel to the local +z- 



axis, we assign +1; if it is anti-parallel, we assign -1. 
In this way, there are spins that are surrounded by the 
equal number of the +1 and the -1. If those spins flip, 
the energy from the first term in Eq.(fT]) does not change, 
but this process is encouraged by the second term. Ear- 
lier numerical calculation confirms that eventually the " 
state" is favored, which occupied +1 and the -1 with the 
equal probability 0,0). Then, this state with the thermal 
fluctuation naturally resolves the inconsistency between 
the experiments. The Edwards-Anderson order parame- 
ter of this state is 0.66770. It contains both the freezing 
and fluctuating sites. Below the freezing temperature 7/, 
the freezing sites contribute to the static internal mag- 
netic field observed in the Ga-NMR experiment, and the 
fluctuating sites remain fluctuating down to zero temper- 
ature to contribute to the reduced moment found in the 
neutron scattering. We note that although the spin vari- 
able on the " + 1" and the "-1" sites may change with the 
coordinate transformation, the one on the "0" site does 
not. The Edward- Anderson order parameter is invariant 
with respect to our coordinate transformation. 

At the mean- field level, as discussed in ref. [J] , the state 
in Fig.([T]) can be stabilized by a ferromagnetic near- 
est neighbor exchange and an anti-ferromagnetic third- 
nearest neighbor exchange. Photoemission spectroscopy 
also supports the enhanced contribution from the third- 
nearest- neighbor sites t 9]. Moreover, in the high tempera- 
ture, the superexchange in this material is believed to be 
isotropic. These complexity concerning the microscopic 
details is not included in this paper. We remark that our 
phenomenological model is built on top of the existence 
of the abnormal spin correlation below the Curie- Weiss 
temperature. The superexchange in our model is anti- 
ferromagnetic in the transformed coordinate. Most im- 
portantly, the existence of the critical phase, as will be 
shown later, and the validity of the current model to de- 
scribe the critical behavior help avoid those complexity. 
In other words, Eq.fl]) might be the fix-point Hamilto- 
nian of the fundamental microscopic model. 

Now, let us focus on the spin-freezing transition ob- 
served at Tf. Since the poly-crystalline sample was used 
in the magnetic susceptibility measurement, we apply the 
quantum Monte Carlo technique to calculate the aver- 
aged susceptibility x defined by l/3(xxx + Xyy + Xzz), 
where Xij = dmi/dhj, and mi is the magnetization per 
site and hj is the external magnetic field. In the cal- 
culation, the 10 6 Monte Carlo Steps with the average 
over 64 ensembles is used. In addition, the cluster algo- 
rithm is applied along the imaginary time direction. In 
Fig-©, the result of the inverse susceptibility x" 1 with 
,/ = 66K (in the temperature unit) and K = 0.5J is 
presented. It is compared qualitatively well with the ex- 
periment and shows weak size-dependence because of its 
average nature. Both of them have the spoon-like shape 
that contains the dip, but the position of the dips happen 
at different temperature from the theoretical result and 
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FIG. 2: (Color online)The system size is 10 x 10 with 100 
spins in the calculation. The theoretical result has the same 
structure as the experiments. At the dip of the x~\ marked 
Tf, the spin- freezing transition occurs found in the Ga-NMR 
experiments. 



the experimental one. We mark the temperature of the 
dip T c in the theoretical result and Tf in the experimental 
one. 
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in the L x x M geometry with the periodic boundary con- 
dition, where L x is the length along the ai direction and 
M is the one along the (1,1,0) direction, and L x ~ 4M is 
taken. In this case, the correlation length £(T, M) scales 
with M. We define 



£(T) = lim 



c(t,m), 



(3) 



and it is shown in Fig.©. We find that the spin correla- 
tion starts to develop at 1.4 J and it diverges at T c ~ 0.3 J. 

At T c , the BKT-type transition occurs. Above T c , for 
example T=0.4J, £(T, M) saturates exponentially shown 
in the Fig.®. Below T c , for example T=0.2J, £(T, M) is 
linear to M. This linearity has only two possibilities: 1) 
the correlation length is so large that our system sizes are 
too small to reach the saturation. 2) it is in the critical 
phase. Because of the conformal invariance of the critical 
phase in the two dimensions, the spin correlation decays 
exponentially in the torus geometry. In this case, £(T, M) 
is linear to M. Furthermore, the scaling dimension A(T) 
can be obtained by 
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FIG. 3: (Color online)The result of £(T). The temperature 
is in the unit of J, and the correlation length is in the unit 
of lattice constant a. At the high temperature side of T c (~ 
0.3 J), the spin correlation is exponentially decayed, and the 
correlation length diverges at T c . Below T c , the phase has the 
power-law spin correlation. The green line is the theoretical 
fit of the 2D XY universality class. 



In the Ga-NMR experiments, the spin-freezing transi- 
tion is observed at the dip of the inverse susceptibility, 
and an unusual spin correlation starts to develop at the 
Curie- Weiss temperature 80K. To understand this, we 
compute the spin-spin correlation along the ai direction 



where A is defined by 



2vrA(T) 



-M 



< a*(0)<r*(r) >~ i 



(4) 



(5) 



The first possibility can be ruled out as the following. 
In Fig.©, we plot the initial slope of £(T, M) defined 
by 9£(T, M)/dM | jvf=o- It shows a clear phase transi- 
tion at T—T c . If the first possibility were true, the ini- 
tial slope would have been monotonic and grew expo- 
nentially as temperature decreases. However, the slope 
has the discontinuity at T c , indicating a phase transition. 
Moreover, it is not a second-order phase transition, be- 
cause both the magnetic susceptibility and the specific 
heat are smooth functions at T c . If it is a classical Ising 
transition, £(T) should be symmetric with respect to T c 
as \T — T c \ << 1. Here, although the initial slopes at 
T=0.2J and T=0.4J are the same within the error bar, 
their asymptotic behaviors are entirely different as shown 
in Fig.®. Therefore, the phase below T c should be criti- 
cal with the power-law spin correlation. Additionally, the 
scaling dimension A(T) has the monotonic temperature 
dependence, which is the typical behavior of the scaling 
dimension in the critical phase. Due to the restriction of 
the technique, we are not able to compute the free energy 
to find the central charge in the critical phase. Classifica- 
tion by the conformal field theory may be an interesting 
directions for the future research. 
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FIG. 4: (Color online)£(T, M) at T=0.2J and 0.4J. The ver- 
tical axis is the £(T, M), and M is in the unit of the lattice 
constant a. At T=0.4J, the correlation length saturates ex- 
ponentially. At T=0.2J, the linear scaling implies the critical 
phase as explained in the text. Two lines are the functional 
fit to the data 
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FIG. 5: (Color online)The temperature is in the unit of J. 
The vertical axis is the initial slope of £(T, M) defined in the 
text. If there is no phase transition, the initial slope should be 
exponentially and monotonically increasing. However, there 
is a discontinuity at T = 0.3 J, which indicates a phase tran- 
sition. 



The existence of the BKT transition in the quan- 
tum Ising model on the triangular lattice was previously 
pointed out by R. Moessner, et al.Q. Here we summa- 
rize their argument and further construct the topological 
object in this model. The quantum 2D model of Eq.|T|) 
at finite temperature can be mapped to a 3D classical 
Ising model with the ferromagnetic exchange along the 
imaginary time direction with finite dimension. Using 
the Landau-Ginzburg- Wilson analysis, the 3D model can 
be mapped to an XY model with a sixth-order symmet- 
ric breaking term which has the sixfold clock symmetry. 
At zero temperature, there is a quantum phase transi- 
tion described by the 3D XY universality class at finite 



K, because the clock term is dangerously irrelevant in 
3D. However, at finite temperature and in the thermo- 
dynamic limit, the 3D model crosses over to the 2D model 
and it results in two transitions: one is BKT transition at 
higher temperature and the other at lower temperature 
corresponds to the sixfold clock symmetry breaking, for 
example, to the state in Fig.fJTJs). Therefore, the BKT 
transition in NiGa2S4 actually belongs to the 2D XY uni- 
versality class! In Fig.f3]), we fit the correlation result 
well with the 2D XY model [13]. These two transitions 
are both seen in the Ga-NMR experiment @. The spin- 
freezing transition at 10K is the BKT transition, and the 
transition at 2K where spins completely freeze up corre- 
sponds to the second transition. 

It is not coincident that the phase transition occurs at 
the dip of the Once the power-law spin correlation 

survives in the phase rather than a critical point, the 
response to the external field could be weaker. Due to the 
slow spin dynamics, the magnetic susceptibility begins to 
drop in the critical phase. 




FIG. 6: (Color online) The temperature dependence of the 
specific heat capacity. The temperature is in Kelvin. 

In Fig.©, we show the calculation of the specific heat 
capacity. It illustrates the double-peak structure and the 
reason is similar to the one given in the Ref. ll| . Al- 
though the valley between the peaks is not as deep as the 
experimental one shown in Ref. [3] • The magnetic specific 
heat in Ref. 0] is obtained by subtracting the specific heat 
of ZnIn2S4, which is non-magnetic and iso-structural to 
NiGa2S4, from the one of the NiGa2S4. We remark that 
there is 29.68% difference in the total atomic mass be- 
tween these two compounds. How reliable their magnetic 
specific result is suspicious to us. However, their low- 
temperature result may be correct. Here we also find 
no evidence of the existence of the energy gap, which is 
consistent with their experiment. 

In summary, we have shown that the novel spin- 
freezing transition seen in the Ga-NMR experiment on 
the NiGa2S4 compound is the BKT- type transition which 
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belongs to the 2D XY universality class. The divergence 
of the spin correlation leads to the divergence of the 1/Ti 
and 1/T 2 . Below Tf, the power-law spin correlation de- 
velops in the state in the Fig.(l). The truly long-ranged 
spin correlation happens at the zero temperature, dubbed 
by the "order by disorder" @, 0] . Through our analysis, 
NiGa2S4 should be removed from the candidate list for 
the spin-liquid ground state. Finally, the long-sought 
BKT transition in the quantum spin system is unexpect- 
edly found. The new phase accompanying with the novel 
phase transition will refresh our understanding of the spin 
dynamics in the critical phase. 
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